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INSTABILITY OF SLENDER THIN-WALLED BOOMS DUE TO 


THERMALLY INDUCED BENDING MOMENTS 
By Vernon K. Merrick 
Ames Research Center 


SUMMARY 


A linearized dynamic stability analysis of slender booms illuminated by 
thermal radiation is developed. It is assumed that the booms have thin- 
walled, circular cross sections. In addition the cross sections are assumed 
to be thermally (but not necessarily structurally) seamless and thermal 
torques about the boom axis are assumed to be negligible. The analysis is 
therefore concerned with the influence of thermal bending moments on boom 
stability. 

It is shown that in the absence of structural damping, any initially 
straight boom, with either a seamless or structurally open cross section, is 
unstable. Booms with an initial curvature may exhibit either oscillatory or 
nonoscillatory instabilities, depending on the magnitude and direction of the 
initial curvature. The worst case oscillatory instability has a frequency 
66 percent of that of the natural flexural frequency and diverges as if it had 
a damping factor of -0.354. The instability can be suppressed by flexural 
damping but, in general, the amount required is greater than can be provided 
by the structure of existing boom types. Torsional damping is relatively 
ineffective and there exists a class of booms which cannot be stabilized by 
torsional damping. 

Initial boom curvature can have a significant effect. If the initial 
curvature is such that the boom, in the equilibrium state, is bent toward the 
radiation source, it is stable; otherwise it is unstable. In general, the 
stability of a boom is only slightly influenced by whether or not the boom 
cross section is restrained from warping at the tip. 

A damper, in the form of a closed vessel, rigidly attached to the boom 
tip and containing a ball free to move through a viscous fluid is an effective 
stabilizer for a large class of booms important in space applications. 


INTRODUCTION 


Consider an idealized situation in which a slender, thin-walled boom, 
clamped rigidly at one end, is placed in a force-free environment. Suppose 
this boom is suddenly illuminated by uniform planar radiation such as, for 
example, solar radiation. Some of the radiant energy will be absorbed by the 
boom and will set up internal thermal gradients. These, in turn, will produce 



internal stresses causing the boom to bend and twist. In the process of bend- 
ing and twisting, elements of the boom will change their orientation relative 
to the direction of the incident radiation. Thus, the thermal stresses and, 
therefore, the effective thermally induced bending moments and torques will 
change as the boom deforms. In concept, at least, the boom will ultimately 
adopt some static equilibrium shape. The calculation of static equilibrium 
shapes for booms of open cross section is treated thoroughly in reference 1. 
However, neither static nor dynamic stability of the boom is considered. There 
is, therefore, no guarantee that all or even any of these static equilibrium 
shapes will occur in practice. 

The first indication that all static equilibrium shapes are not stable 
was both sudden and dramatic. The 0G0 IV satellite which has a 60-foot- long 
boom of open cross section began to exhibit high frequency attitude oscilla- 
tions whenever it was exposed to sunlight. It was found by the Stabilization 
and Control Branch at Goddard Space Flight Center that this behavior could be 
explained if it was assumed that the boom was oscillating at its fundamental 
bending frequency with an amplitude at the tip of about 14 feet. This infor- 
mation became available at a time when many gravity stabilized satellites were 
exhibiting apparently random anomalous attitude behavior. Since long booms of 
open cross section are common features of these satellites, it was conjectured 
that boom instability could be the source of the difficulties. High frequency 
oscillations recently observed in the magnetometer readings of one of these 
satellites tend to support this conjecture. The final evidence supporting the 
possibility of boom instability was provided by the Structural Dynamics Branch 
at Ames Research Center where thermal oscillations were produced in a rela- 
tively short boom. This demonstration provided incontrovertible proof of the 
possibility of thermally induced dynamic instability in booms of open cross 
section. 

In retrospect, a hint that the stability of booms may be a significant 
question is given in reference 1, which points out, clearly, that the bending 
and twisting modes are coupled. The existence of coupled structural modes 
together with bending moments and torques that are functions of structural 
deformation are classical elements in the production of instabilities. 

Numerous examples can be found in the field of aeroelasticity . 

There are several approaches to the theoretical study of thermally induced 
boom instabilities . An extension of the comprehensive approach of reference 1 
to include dynamic effects should produce details of the actual boom motions 
that will be encountered. An important contribution along these lines has been 
given by Frisch in reference 2. However, with this type of approach it is 
often difficult to assess the relative importance of the various boom system 
parameters. Therefore, in addition to a precise analysis, there is a need for 
a more traditional engineering approach in which the model of the phenomenon 
is drastically simplified to the point where only the most important factors 
are retained. With such a theory there is hope of obtaining insight into the 
basic causes of boom instability and into possible means of preventing it. 
Analysis of possible thermal instability mechanisms have already been made by 
Yu (ref . 3) and Beam (ref. 4) . The mechanisms in both cases were relatively 
simple although they differed widely in concept, Yu analyzed a situation in 
which changes in thermal bending moments, due solely to boom flexure, provide 
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the mechanism for converting heat energy to mechanical energy. The type of 
boom cross section in this case is not important. One result of this theory 
is that the greatest instabilities occur when the axis of the boom lies along 
the direction of the radiation. In contrast. Beam presents an analysis, sup- 
ported by experiment, of the influence of thermal torques on stability. An 
essential feature of the type of boom treated is that it has an open nonover- 
lapping cross section. The greatest influence of thermal torques must, of 
course, occur when the boom axis is perpendicular to the direction of radia- 
tion. Beam discovered that booms could exhibit both instability and improved 
damping, depending on the orientation of the boom split relative to the 
direction of radiation. 

The purpose of this report is to analyze yet another mechanism for the 
production of thermal instability. Consider a boom which, in its unheated 
state, has its axis normal to the direction of radiation. The static equilib- 
rium shape will generally be such that the boom is bent away from the radia- 
tion source. A bent boom has the well-known structural property of coupled 
flexure and torsion. Thus a force applied to the end of the boom in a direc- 
tion perpendicular to the plane containing the boom axis will produce not only 
bending but also twisting. This effect is particularly marked for a boom of 
low torsional rigidity. The boom twist produced by the force causes thermal 
bending moments that tend to restore the plane of bending back to its original 
direction. The analysis presented here seeks to determine how boom stability 
is influenced by bending moments due to twist. It is clear that the basic 
physical mechanism described above has many similarities to that underlying 
the in-flight dynamic behavior of an airplane wing. Once this analogy is 
recognized it becomes easier to understand how a thermally bent boom can 
exhibit instabilities. Since in a practical boom system employing an open 
cross section thermal torques and thermal bending moments due to twist must 
coexist, the analysis of this report may be regarded as complementary to that 
of Beam in that any general theory of boom stability must contain essential 
elements from both theories. 

The analytical approach adopted here is based on the linearized analysis 
of what is probably the simplest conceivable boom system. One major simplify- 
ing characteristic of this system is derived from the assumption that the 
cross section provides a continuous conduction heat path. The cross section 
is therefore thermally seamless (although not necessarily structurally seam- 
less), which greatly simplifies the calculation of the thermal bending 
moments. A second major simplifying characteristic is derived from the 
assumption that the effective mass of the boom may be considered to be located 
at the free end. This assumption ensures the existence of a single mode of 
bending and a single mode of twist. A third major simplification is that 
thermal torques are negligible. This assumption is strictly true only if the 
cross section is structurally seamless. This boom system is not particularly 
representative of the open overlapped cross section booms used in space 
applications. On the other hand, the behavior of the idealized boom system 
can provide basic information relevant to an understanding of the observed 
behavior of actual boom systems. In particular, the analysis provides an 
immediate answer to the basic question of whether or not thermal bending 
moment instabilities are possible with a boom having a structurally seamless 
cross section. 



The analysis is extended beyond that which would be required to analyze 
the stability of a simple, initially straight boom. The influence of small 
amounts of initial boom curvature (such as might be introduced during manu- 
facture) is investigated, along with the possibility or providing stability by 
means of a simple, passive, ball damper located at the tip of the boom. 

DERIVATION OF THE STABILITY POLYNOMIAL 


Four steps are required to develop the stability polynomial. The first 
step is to describe geometrically the static equilibrium shape of the boom in 
terms of its assumed initial, unheated shape and the bending due to thermal 
radiation. This static equilibrium shape is described in terms of a radius of 
curvature and the. angle between the final plane of bending and some fixed ref- 
erence. These variables appear in the stability polynomial. The second step 
is to establish a relationship between the thermal bending moments and the 
angle of twist of the boom. The details of this step are given in appendix A. 
The third step is to determine the structural properties of the boom. These 
take the form of load-deflection relationships. The loads are assumed to be 
forces acting at the boom tip in a direction normal to the boom axis, a torque 
acting along the boom axis at the tip, and an arbitrary distribution of bend- 
ing moments along the boom axis. The appropriate equations are developed in 
appendix B. The final step is to identify the loads acting on the boom in 



Figure 1.- Boom system (primary coordinates and 
geometry of deformation) . 


terms of force of inertia of the tip 
mass and the thermal bending moments. 
Four linear differential equations 
result whose characteristic polynomial 
defines the stability of the boom 
system. 


Principal Assumptions 

The boom system to be analyzed and 
the primary coordinate reference system 
is shown in figure 1. A precise 
definition of this coordinate system is 
given in the Nomenclature. 

The following principal assump- 
tions are made about the characteristics 
of the boom. 

1. In the absence of any radia- 
tion the centroidal axis of the boom 
has a constant radius of curvature R Q 
along the length of the boom. Further- 
more, the curvature vector at any point 
on the boom axis lies in a plane whose 
normal makes an angle 6 to the y Q r 
axis (see fig. 1). 
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2. The initial radius of curvature R 0 and the thermally induced 
radius of curvature R s are both large compared with the boom length L. 

3. The boom cross section is circular, thin-walled, thermally seamless, 
and constant along the length. In addition, the principal second moments of 
area are equal. 

4. Warping of the boom cross section is restrained at the root. 

5. The entire mass of the boom system is located at the free end. In 
practice this means that the tip mass of the boom is large compared with the 
mass of the boom itself. 

6. Boom bending does not substantially change the orientation of the 
boom relative to the direction of radiation and, therefore, has no influence 
on thermal bending moments. This assumption places some restrictions on the 
magnitude of the angle 3 between the y Q ! axis and the direction of the 
radi ation. 

7. The centroidal and shear axes of the boom coincide. 

8. Thermal torques are negligible. 

Several other special assumptions of a less fundamental nature are 
introduced and defined as they are required in the .analysis. 


Geometry of the Equilibrium State 

The static equilibrium shape of the boom results from the combined 
effects of its initial shape and a thermally induced bending moment. It 
follows from equations (A29) that the thermally induced bending moment is con- 
stant along the boom axis and acts about the z Q ? axis. Simple bending 
theory shows that the thermally induced equilibrium radius of curvature R s 
is constant along the length and the corresponding contribution to the total 
curvature vector lies in the x^y^ plane. It follows that, to a first 
approximation, the static equilibrium shape resulting from the combined 
effects of initial and thermally induced curvatures is planar. This is 
because for all points along the boom the triangles OPQ, shown in figure 1, 
are similar. In fact, at any point x along the boom OP = x 2 /2R Q , 

PQ » x 2 /2R s , OQ a x 2 / 2R, where R is the radius of curvature of the static 
equilibrium shape, and R Q is the radius of curvature of the initial shape. 
The geometry shown in figure 1 then yields the following expressions from 
which the radius of curvature R and the angle a between the normal to the 
plane of the boom and the y Q 1 axis may be obtained. 

_ J_ 

R 2 “ R 0 2 + 


Rc 


2 sin 6 

R 0 R S 


CD 


i 
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sin a 
R 


sin 6 


( 2 ) 



R 
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Thermal Bending Moments and Boom Deflection Equations 


Suppose now that the boom system 
forces Fy(t) and F z (t) and a torque 



is deflected from its equilibrium by 
T(t) (see fig. 2). By assumption (6) 
the translational deflections y and z 
cause no changes in the thermal bending 
moments . On the other hand it is shown 
in appendix A that a twist deflection ip 
does produce changes in the thermal 
bending moments . It is shown in appen- 
dix A that for small angles of twist, 
the equilibrium thermal bending moment 
M ,,(») remains constant in magnitude. 

However, if it is assumed that the 
direction of the M ,, (°°) moment rotates 

with the cross section, then an addi- 
tional thermal bending moment 
AM , (x,t) whose magnitude is given by 

equation (A37) acts in the direction of 
the y' axis. Thus, 

d AMy t (x,t) 

^ + A AMy. (x,t) = AM z „ (°°)iKx,t) 

(3) 

and from equations (A29) and (A15) 




J s a s e c Eirr 2 cos 3 
2p cA 


(4) 


4gcT 0 3 + K 
hpc pcr2 


(5) 


where B is the angle between the direction of radiation and the negative 
y Q ' axis (see fig. 1). (Note: since thermal bending moments are associated 
only with the first thermal mode, the subscript may be deleted from Aj.) 

Equation (3) may be written in the following form 


dn(x,t) , 
dt 


An (x,t) 


Ai|i(x,t) 


( 6 ) 
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where 


I 


n(x,t) 


AM y , (x,t) 

M Z „C°°) 


(7) 


Since (°°) remains constant in magnitude and its direction rotates with the 

cross section, it follows that its entire effect is equivalent to a permanent 
bend in the boom. The radius of curvature of the bend, as given by the 
Bemoulli-Euler bending theory, is 



m z „M 

El 


( 8 ) 


With M !f (°°) replaced by this permanent bend, the only thermal bending moment 
acting on the boom is AM f (x,t). Since by assumption (8) thermal torques are 

assumed to be negligible, the load system acting on the boom is as shown in 
figure 2. The deflections of the boom may be obtained directly from the 
results of appendix B. It follows from figure 2 that 


AM^(x,t) = -AMy f (x, t) cos a = - -g— n(x,t)cos a 

El 

AM n (x,t) = AM f (x,t)sin a = n(x,t)sin a 
y k 


(9) 


Incorporating the above relationships into equations (B34), (B35), and (B22) 
gives the following expressions for y(L,t), z(L,t), and ^(L,t): 


yCb,t) 


Fy(t) 

3EI 




N(L,t)cos a 


R 
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( 10 ) 


z(L,t) = FJt) 


{l 3 l 5 

rfa(kL) 
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1 , T (t) L 3 
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1 1 

\ 3EI 5R 2 

4C 
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3R 

2C 

" ElJ 


N (L,t)sin a 
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(ID 


iKL,t) 


F z Ct)L 3 rfjfkL) 
3R 2C 


1 H T(t)Lf 2 (kL) N(L,t)sin a 
El J C RR S 


where 


N(L,t) 



n (x,t)sin 



dx 


( 12 ) 


(13) 
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and graphs of the functions fi(kL), 
f 2 (kL), and f 3 (kL) are given in 
figure 3. 

Equation (6) is now multiplied by 
R sin[(L - x)/R] and integrated with 
respect to x over the interval 
0 <_ x <. L. The result is 

^ dt + AN(L,t) = AR^|* iKx,t)sin^^^dx 

(14) 


Figure 3.- Variation of f x (kL), f 2 (kL), £ 3 (ki) , It has been shown in appendix C 

and f 4 (kL). that over any interval of time starting 

at zero the maximum angle of twist due to thermal bending moments alone is 
less than the maximum total angle of twist multiplied by (x 2 /2RR s )sin a. 

Since, by assumption (2), (x 2 /2RR s ) << 1, it follows that angle of twist due 
to thermal bending moments may be neglected without much loss of accuracy. 

Thus the last term of equation (12) may be deleted along with the last term of 
equation (B21). This simplification of equation (B21) permits the integral 
expression in equation (14) to be evaluated explicitly. Thus, the integral 
can be written in the form 


R f ip(x,t)sin 

•'O 




+ 


4* b (x, t ) ] sin 



(15) 


where i^(x,t) and 4;g(x,t) are defined in appendix B. Therefore by 
equations (B31) and (B32) , equation (15) becomes 


/' 


i^(x,t)sin 
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dx = 


F z (t)L 5 f q 


20 RC 


CkL) T (t) L 
+ — 


‘fi (kL) 


6C 


+ R 


/' 


tJ> B (x,t)sin 


(t) 


dx 


(16) 


Substituting for i|/g(x,t) from equation (B21) into equation (16) gives the 
following final form of equation (14) 


dN (L, t) , 
dt 


AN (L, t) 


Equations (10), (11), (12), 
deflection equations. 


AF z (t) L5 rf 3 ( kL ) i "I AT(t)L 3 f! (kL) 

10R [ 2C 3EIJ + 6C 

and (17) are the final complete set of boom 


(17) 
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The Characteristic Stability Polynomial 


If the boom system has no external forces acting on it, the forces F (t), 
F z (t) and torque T(t) are due entirely to mass -accelerations and internal 7 
structural damping and may be represented as follows: 


F ft i - n, d 2 y(L,t) dy(L,t) 

F y(t) - -ittj dt2 - D y dt 

(18) 


F ftl - m d2z t L » t 5 n dz(L,t) 

F z Ct) ~ " m T dt 2 °z dt 

(19) 


xm _ T d 2 iKL,t) n di|>(L,t) 

UtJ dt 2 ^ dt 

(20) 

where 



mrp 

tip mass 


s 

tip inertia about the boom axis 


D y ,D z ,c V 

effective structural damping for translation along b(L) and 
and rotation about t(L), respectively 

n(L) 


Equations (10), (11), (12), (17), (18), (19), and (20) form the complete set 

of linearized dynamical equations from which F (t) , F z (t), and T(t) may be 

eliminated to yield four equations in the variables y(L,t), z(L,t), ip(L,t), 
and N(L,t). These equations describe the motion of the boom system under the 
sole influence of the thermoelastic effects due to thermal radiation. The 
characteristic polynomial is given below in determinant form where the symbol 
"s" is used to denote the Laplace transform variable. The columns of the 
determinant correspond to the variables y(L,t), z(L,t), ip(L,t), and N(L,t), 
respectively. The first two rows correspond to the y and z translational 
equations of motion of the tip, the third row to the rotational equation of 

motion of the tip about the boom axis and the fourth row to the dynamic rela- 

tionship between thermal bending moments and angle of twist. 


(m T s 2 + U y s)(l - ^2*) 


(l 3 L 5 

Cm T s 2 + D z s) 3 e T + ^2 


[~ f 3 (kL) _ _j_ll 

l 4C 3ElJj 

<-r- ! * & 

<t*’ * “.‘Jp# 1 - m] iro 


4- 1 


„ s 2 + d of flCkL) , JM Lj 
U i^ s % SJ 1_ 2C El J 3R 




(I*s 2 + V 5 


L 3 f j (kL) A 
6C 


cos a 
R c 


= 0 


S + A 


( 21 ) 


Mi 



One solution of equation (21) is 


Cm T s* * D y s)(l - 1 = 0 (22) 

This is a damped mode of oscillation for any positive values of D^. It 
cannot exhibit instability and is therefore not considered further. The 
stability of the remaining modes of motion is given by the remaining solutions 
of equation (21) : 


. ( L 3 l 5 

("h 8 ‘hsJjjITV 


(V 2 + V)[- 

l(kL) 11 
2C El J 

L 3 

3R 

sin a 
R s 



(TS 2 * D z s)[ ) 


(Hs 2 + V) 

Lf 2 (kL) 

c 

1 

0 

= 0 

(23) 

(V 2 + »,»)[ 3 2C ) 

l ] l 5 a 
3EI J 10R 

CV 2 * V) 

L 3 f!(kL)A 

6C 


s + A 




ANALYSIS OF THE STABILITY POLYNOMIAL 


In this section the stability polynomial, given by equation (23), is 
analyzed to obtain the maximum degree of instability and conditions for sta- 
bility. Two types of boom cross section are considered. The first is the 
structurally seamless cross section which has the important characteristic 
that the torsional and flexural rigidities are of the same order of magnitude. 
The ratio of the two is given by the expression 


C 1 

El ~ 1 + v 


(24) 


where v is Poisson ! s ratio. For most materials Poisson’s ratio is close to 
1/3 and, when specific numerical results are given, this is the value that 
will be assumed. The second cross section considered is the open overlapped 
type used in most space applications. The important structural characteristic 
of this cross section is that the torsional stiffness is much smaller than the 
bending stiffness. This may be expressed in the form 

Wt <k 1 (25) 


System Stability in the Absence of Structural Damping 
When = D z = 0, equation (23) can be written in the form 


10 



1 


1 


where 


mjs 2 A + 
m T s 2 D 
nvp s 2 B A 


V 2d 

I^s 2 F + 1 0 

I^s 2 GA s + A 


0 


A A 


L 3 L 5 ff 3 (kL) 

3EI 5R 2 L 4C 



B 4 
D k 


f 3 (kL) 

2C 

'fi(kL) 

2C 


1 1 L 5 sin a 
3EIJ ~10RR S 


1 1 L 3 
EIJ 3R 


F = j=r f 2 (kL) 

^ A L 3 fi(kL)sin a 

^ — r n r> 


Equation (26) , in its expanded form, is 


s 5 m T I l(j (AF - D 2 ) + s 4 m T I^A[D(G - D) + F (A - B)] 

+ s 3 (m T A + I^F) + s 2 A[m T (A - B) + I^F] + s + A = 0 

It is shown in reference 5 that any fifth-order polynomial 

a Q s 5 + a^s 4 + a2S 3 + a3S 2 + a4S + as = 0 

has stable roots if and only if all the coefficients are the same sign 
following two inequalities are also satisfied simultaneously 

a! a 2 - a 0 a 3 > 0 ' 

(a 1 a 2 — a 0 a 3 ) (a^a^ — a 2 a 3 ) - (a^a^. — aga 3 ) 2 > 0 ( 

It follows from equation (32) that 


(26) 

(27) 

(28) 

(29) 

(30) 

(31) 

(32) 

and the 

(33) 
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3-33-4 ~ a 2<*5 “ -m^BX 


Therefore when B > 0 


a 3 a 4 - a 2 a 5 < 0 (34) 

However, if condition (34) holds then conditions (33) cannot hold simulta- 
neously. Therefore the system is unstable. In particular, for initially 
straight booms of either seamless or open cross section B > 0. These 
initially straight booms are therefore always unstable in the absence of 
structural damping . 


Conditions Under Which 1^ May Be Neglected 

Consider the case of a boom not under the influence of thermal radiation 
but still bent into a circular arc of radius R. The coupled flexural- 
torsional dynamical behavior of this boom would be described by the transla- 
tional equation of motion in the z direction and the rotational equation of 
motion about the boom axis. Since heating is absent, the variable N and the 
thermal bending moment equation need not be included. The frequencies of the 
natural undamped modes of vibration (D z = = 0) are therefore given by the 

first two rows and columns of the determinant in equation (26); 


m^As 2 + 1 
m T Ds 2 


I^Ds 2 
I^Fs 2 + 1 


(35) 


Equation (35) may be expanded into the form 


m-pI^AF - D 2 ) s 4 + (I^F + m T A)s 2 + 1 = 0 


(36) 


If 1^ is to have an insignificant effect on the frequency of the mode 
associated with m-p and which, in an unbent boom, would be the flexural mode, 
the polynomial m-pAs 2 + 1 must be an approximate factor of equation (35). It 
can be shown, by dividing equation (36) by m-pAs 2 + 1, that this is true if 
and only if 



1 


A good working rule, therefore, might be that I. has very little influence 
on the flexural frequency if ^ 


I 




(37) 
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If is the radius of gyration of the tip mass about the axis of the boom, 

then inequality (37) can be written in the form 


P 



A 

D 


(38) 


It will be demonstrated, in the sequel, that the influence of thermal radia- 
tion is to produce an unstable mode with a frequency close to that given by 
m T As 2 + 1 = 0. It is also clear that reasonable amounts of structural damping 

will not have a great influence on the frequency of this unstable mode. There- 
fore inequality (38), which in dicates when 1^ does not influence the natural 
structural mode of frequency /l/m^A, should also give a rough indication of 
when 1^ does not influence the thermally induced instability of a boom 
possessing some structural damping. 

For a boom with a structurally seamless cross section k = °°, so that 
fi (kL) = f 2 (kL) = f 3 (kL) = 1, and the expressions for A and D given by 
equations (27) and (29) become (with v = 1/3) 


A = 


L 3 

3EI 


D 


L 3 

9REI 


Inequality (38) then becomes 


. 3R 

P , < — — 

* " /lo 


* R 


(39) 


Most booms have a value of R greater than 100 m so that it is clear from 
inequality (39) that, for most practical boom systems employing structurally 


seamless cross sections, the tip inertia 
boom dynamics. 


1^ will have a negligible effect on 


For a boom with an open cross section, A and D become 


D = 


L 5 £ 3 (kL) 


3EI 20R Z C 

L 3 fx (kL) 

6RC 


Inequality (38) then becomes 
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(40) 


L 3 . L 5 f 3 (kL) 
< 1 3EI 20R 2 C 

* ~~ /l0 L 3 £i (kL) 

6RC 


The requirements on given by inequality (40) are most stringent when 

C -*■ 0. Also the ratio f 3 (kL)/f 3 (kL) is close to unity for most booms. Thus 
1^ cannot influence the dynamics of booms with open cross section (C/EI << 1) 
if 


P 


i|> - 


3 \J_ 

10/10 R 


L 2 

10R 


(41) 


It can be seen by comparing inequalities (39) and (41) that the limiting value 
of for a boom with an open cross section is much less than that for a 

boom with a seamless cross section. In fact, the ratio of the two is 
(1/10) (L/R) 2 . Even though the limiting value of indicated by inequal- 

ity (41) can be quite low, it has never been exceeded in boom systems used for 
space applications. Consequently, for the remainder of the analytical treat- 
ment given in this report it is assumed that 1^ = 0. It is conceivable, 
however, that future boom systems may employ tip masses with rotational 
inertias large enough to influence the boom dynamics. In these cases the 
complete polynomial, given by equation (23), must be used in evaluating the 
stability. 


The Influence of Flexural Damping D z 

Some of the implications of equation (23) can be best understood by a 
study of the system with no rotational inertia or rotational damping. With 
1^ = = 0 equation (23) reduces to 


m T As 3 + [D Z A + m T X(A - B)]s 2 + [1 + D Z X(A - B)]s + X = 0 (42) 

For a seamless cross section, equation (24) and L 2 /R 2 << 1, permit equa- 
tions (27) and (28) to be simplified: 


A = 


L 3 

3EI 


(43) 


(1 + 3v)L 5 sin a 
60EIRR S 


f 3 (kL) = 1 


For a split cross section using inequality (25), the corresponding 
expressions are 


(44) 
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+ 


A = 


L 3 

3EI 


L 5 f 3 (kL) 

20CR 2 


(45) 


L 5 f 3 (kL)sin a 
2 OCRR 3 


(46) 


From the points 
results , a more 
in terms of the 


of view of both analysis and presentation of computational 
convenient form of equation (42) is obtained by expressing it 
following set of dimensionless parameters: 


X* k x/m T (A - B) (47) 

s* ^ s/m T (A - B) (48) 

c 2 4 /^(A - B) (49) 

X i j (50) 


It should be not ed that c is the flexural damping factor based on a circu- 
lar frequency /m^(A - B) . Normally the flexural damping factor t z f is 

based on the natural undamped circular frequency in flexure /m T L3/3EI . The 
conversion from one to the other is therefore 




*z/ 3EI (A - B) 


The dimensionless form of equation (42) corresponding to the dimensionless 
parameters defined above is 

* 3 2c 

x 5 * x + ( j ~ + A*)s* 2 + (1 + 2X*C z )s* + A* = 0 (51) 

It is shown in reference 5 that any cubic equation 


a Q s 3 + ais 2 + a 2 s + a 3 = 0 


has stable roots if and only if all the coefficients have the same sign and 

aia 2 > a 0 a 3 (52) 

Inequality (52) written in terms of the coefficients of equation (51) becomes 


hi 
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(53) 


e 


2 


Z 





[1 + A * 2 (1 - X)] - J > 0 


If no structural damping is present (c z = 0) the stability criteria can be 
satisfied if and only if X < 0. It follows from the definition of X that 
this condition is equivalent to sin a<0 or 0 > a > -tt. Therefore stabil- 
ity or instability will prevail according to whether or not the boom, in 
static equilibrium, and the radiation source are both on the same side of a 
plane passing through the boom root perpendicular to the direction of radia- 
tion. This condition for stability may be interpreted in terms of a relation- 
ship betv/een the direction and magnitude of the initial curvature and the 
thermally induced curvature. Thus, from equation (2), if (sin 6 /R Q ) < - (1/R S ) 
the boom will be stable and if (sin 6/R 0 ) > -(1/R S ) the boom will be unstable. 
In particular, if the boom is initially straight, then (1/R 0 ) = 0, and it is 
always unstable. The instability can occur in two forms, depending on the 
value of X. This can be seen more clearly if equation (51), with ? z = 0, is 
expressed in the standard root locus form 

i + [a* ( i - x)] — = o 

s*(s* + 1 - X) 


in which A*(l - X) is the gain factor. 


The root locus for 0 < X < ] 
be seen from figure 4. Thus when 

Imaginary Imaginary 


t> 


P 




differs from that for 1 < X < <*> as can 
0 < X < 1 the instability is always an 
oscillatory divergence; whereas when 
1 < X < 00 the instability may be a pure 
divergence of the type associated with 
static instability. 


(a) 0 < x < 1 


R eai 1 | Real It follows from equations (43) and 

(44) that the value of X for a boom 
with a seamless cross section can never 
exceed [(1 + 3v)L 2 sin a]/20RR s , which 
(b) 1 < x < co is sma H compared with unity. For such 
a boom, therefore, any instability must 
always be of the oscillatory or flutter 
type. This same statement, however, 

cannot be made for booms of open cross section. In this case it follows from 
equations (45) and (46) that 

3 EI L 2 f 3 (kL) sin a 
20 C RR S "" 


Figure 4.- Root Loci for X > 0. 


X = 


-i 

1 + 


3 EIL 2 f 3 (kL) 
20 CR 2 


(54) 


If 1/R 2 and sin ot/R, given by equations (1) and (2), respectively, are sub- 
stituted into equation (54), there results an expression for X in terms of 
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R 0 and 6. The maximum value of X with respect to R 0 and 6 is given by 


y _ 1 /3EIL 2 f 3 (kL) 

2R S / 20C 

This maximum occurs when 6 = — it/ 2 and 


1 1 / 20C 

R q R s / 3EIL 2 f 3 (kL) 


Clearly, therefore, as EI/C X -*■ <*>. However, from equation (54), if the 

values of R Q and 6 satisfy the inequality 

1 / sin 6 + 1 \ > -2QC 

R o \ R s + R o / 3EIL 2 f 3 (kL) 

then X <_ 1. This condition on R Q and 6 is not too stringent and, in 
particular, holds for initially straight booms. It is anticipated that most 
practical boom systems will satisfy the above condition on R Q and 5, which is 
therefore assumed to hold throughout the remainder of this report. The report 
is therefore concerned only with the oscillatory instabilities associated with 
values of X in the interval 0 < X ^ 1. 

It is of interest to know the greatest degree of instability that can 
occur, as measured by magnitude of the real part of the divergent roots of 
equation (51). The answer to this question depends on which of the two vari- 
ables A*, X is maintained constant. Both cases, namely, variable A* at a 
fixed X and variable X at a fixed A* are treated in appendix D. When 
C z = 0 and X is constant, the real part of the divergent root has a maximum 
value given by 

P* = f (55) 


and this occurs at a value of A* given by 

x* - 2 - x 

A 2(1 - X) 


(56) 


with the corresponding imaginary part of the root given by 


q* = 



X xi 
2 " 16 


(57) 


111 
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These results may be expressed in dimensional form and simplified for the case 
of a boom with a seamless cross section by substituting for A and B from 
equations (43) and (44) and using the condition that X << 1. The results, in 
dimensional form (see notation for definition of p* and q*), are 


L 2 (l + 3v)sin a / 3EI 
80RR^ / m T L 3 


(58) 



Equations (58) and (59) show that since L 2 sin a/RR s << 1 the frequency of 
the divergent mode is approximately equal to the frequency of the flexural 
mode of a straight boom. The instability will therefore appear to an observer 
as if the boom were oscillating in both bending and torsion at its natural 
bending frequency, and diverging as if it had a damping factor equal to 
[-L 2 (l + 3v)sin a]/80RR s . If 1/Ro > 1/R S this quantity has a lower bound of 
[-(1 + 3v)/40] (L 2 /R s 2 ) . With v = 1/3 and L/R s =0.1 this lower bound has a 
value of -0.0005. This is clearly a mild divergence. The basic reason why 
booms with seamless cross section can exhibit only a mild divergence is that, 
because the torsional and bending stiffnesses are of the same order of magni- 
tude, X = B/A is always small. This is not the case when the torsional 
stiffness is low compared with the bending stiffness, as for a boom of open 
cross section. Thus, for example, equation (54) shows that, for an initially 
straight boom (1/R 0 = 0), X 1 as EI/C °°. 

When £ z = 0, it follows from the results of appendix D that, for a 
fixed value of X*, the maximum value of the real part of the divergent root 
is (eq. (D29) ) 


P 2(1 + 2X*) 

and this occurs at a value of X given by (eq. (D30)) 


(60) 


X = 1 


X* 

(i + x*rcr + 2x*) 


(61) 


with the corresponding imaginary part of the root given by (eq. (D32) using 
eqs . (60) and (D28)) 


q* 


/X* (4 + 7X*) 
4(1 + 2X*) 2 


(62) 


18 



Maximum p 


I 



CL 

E 

3 

E 

X 

o 

E 


X 


Graphs of equations (60), (61), and (62) 
are given in figure 5. 

It follows from equations (60), 
(61), and (62) that p* attains its 
maximum value when A* «>. In the 
limit 


lim p* = j (63) 

A*+~ 


lim X = 1 (64) 

A*-* 00 


Figure 5.- Variation of maximum real part of 
unstable root (p*) with X*. 


lim q* = /L (65) 

A*^ OO 


For the case of an initially straight boom (1/R Q = 0) the corresponding 
limiting dimensional values of the unstable roots are 


I /3EI~ 
^ 4 / impL 3 


( 66 ) 


q 


HT_ /5EI" 
v 16 v mjL3 


(67) 


It follows from equations (66) and (67) that 



( 68 ) 



Figure 6.- Variation of real part of unstable 
root (p*) , with X*, for various values 
of X. 


where /p 2 + q 2 is the circular fre- 
quency of the unstable oscillation. To 
an observer, therefore, a boom for which 
X 1 and A* + 00 will appear to oscil- 
late in both bending and torsion with a 
frequency that is about 66 percent of the 
natural bending frequency of a straight 
boom. It will also appear to be diverg- 
ing as if it had a damping factor equal 
to -v/2/4 or -0.354. This is clearly a 
rapid divergence. Furthermore, it 
remains severe down to quite a small 
value of A*. For example, it is shown 
in figure 5 that if A* > 2, then 
p* > 0.2. To complete the survey of the 
instabilities shown by equation (51) the 
variation of p* with A* for various 
values of X is shown in figure 6. This 
graph shows the two types of maxima for 
p* which have been discussed. 
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It is important to note here that the effect of boom tip condition on 
stability is dictated by the function fs(kL) (see eq. (54)). It follows from 
figure 3 that, for values of kL greater than about 2, the stability of a 
boom is independent of whether or not the tip is restrained from warping. 

For a fixed value of X the amount of flexural damping ? z required to 
stabilize the boom for all values of A* may be determined by setting p* = 0 
in equation (D16) . The result is 


?z = J ' J /rrY ( 69 ) 

and occurs at the following value of A* obtained by substituting equa- 
tion (69) into equation (D17), 


A* = - — (70) 

/I - X 

with the corresponding imaginary part of the root given by 

q* 2 = /! - X (71) 


Alternatively, equations (69) and (70) may be obtained by differentiating 
equation (53) with respect to A* and setting 3? Z /3A* = 0. 

For a boom with a seamless cross section equations (69), (70), and (71) 
become, using equations (43) and (44), 


L 2 (l + 3v)sin a 
80RR S 


(72) 


A 


q = 



(73) 


These results compared with equations (58) and (59) show that the structural 
damping required to counter the boom instability is equal to the magnitude of 
the effective negative damping factor corresponding to the maximum instability 
when the boom has no structural damping. The value of £ z given by equa- 
tion (72) is very small and will probably be exceeded by the structural damp- 
ing present in most practical boom systems employing seamless cross-section 
members. It is, therefore, unlikely that instability will ever be observed in 
these types of boom. 

The maximum structural damping required for stability when A* is 
assumed constant occurs when X = 1. This can be shown from equations (D29) 
and (D30) of appendix D, from which it also follows that 
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-1 + 


A 


+ 4A* 


4A' 


( 74 ) 


for p* = 0. The corresponding value for 
the imaginary part of the root is q* = 0. 
In other words the oscillation is M dead 
beat." The graph of equation (74) is 
plotted in figure 7 along with the varia- 
tion of with A* for other values of 
X. Since for a boom with an open cross 
section X can approach unity, it is 
clear that large values of £ z may be 
Figure 7.- Minimum flexural damping U z ) required required. In particular, if A* > 2 

for stability. then £ z > 0.4. The asymptotic value of 

for X = 1 is 0.5. This therefore represents the c z which will stabi- 
lize a boom with a circular thin-wall cross section irrespective of the mate- 
rial it is made from, the initial curvature, or the L/R s ratio. It is 
interesting to note that the £ z required to stabilize the worst case insta- 
bilities is somewhat greater than the negative damping factor of the 
instability when £ = 0. 


The Influence of Torsional Damping 

It has been demonstrated in the last section that it is unlikely that 
instabilities of booms with seamless cross sections will ever be observed. 
The following treatment, to determine the influence of torsional damping, is 
therefore restricted to the more unstable open cross-section boom configura- 
tions characterized by EI/C >> 1. Consider the case wherein D z = 1^ = 0 
and the boom is initially straight (1/R 0 = 0)* The stability polynomial 
(eq. (23)) then reduces to 


2 [9 - 5X 0 f 4 (kL) l *4 ( 1 A *3 

9 SL 1 - X Q J S ° \1 - X 0 ^<F° ) ° 

+ (A 0 * + 2 ^) s q * 2 + (2^A 0 * + l)s Q * + A Q ‘ 
_ D^Lf 2 (kL) /SET 


where 


= 0 




2C 


m^L 3 


3 EI_ ]J_ _ . 

20 C R2 ^3 ( kL ) 


, '3 El L2 

1 + 20 ~C ^3 ( kL ) 


v * _ E 

V 3EI 


(75) 

(76) 

(77) 

(78) 
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(79) 


f 4 (kL) 


[fl (kL)3 2 

f 2 (kL)f 3 (kL) 


Note that X 0 , A Q *, and s Q * are the values of X, A*, and s*, respectively, 
for an initially straight boom. The graph of f 4 (kL) is given in figure 3. 

It is shown in reference 5 that any fourth-order polynomial 

a Q s 4 + ais 3 + a 2 s 2 + a 3 s + a 4 = 0 

has stable roots if and only if all the coefficients are positive and 

a 3 (ai a 2 - ao a 3) _ a 4 a l 2 > 0 (80) 

In this instance this condition reduces to 


8 [9 - SfaklOJXo* 2 ?^ 3 + 4 [9 - 10f 4 (kL)]A o *^ 2 


9X * 

- [10f 4 (kL) - 18A q * 2 ]^ + — — - < 0 


(81) 


\ * 
A o 


space 


-c* 


f 4 CkL) 


(Note Boom is initially straight) 


This inequality defines a stability boundary which divides the X Q , 
into two regions. These are characterized by whether or not the stability 
condition can be satisfied by a real positive value of If the stability 

condition can be satisfied, then the system can be stabilized by torsional 
damping; otherwise it cannot. This stability boundary has been calculated for 

several values of ft+CkL) and is shown in 
figure 8 with the values of along 

the boundary. The cubic equation in 
defining the boundary has two equal real 
positive roots and one real negative root. 
In the left hand region of figure 8 there 
are usually two real positive, but 
unequal, roots and one negative real root. 
The boom system can then be stabilized 



ft only for values of ^ in the closed 


Values of 


for which the 


Figure 8.- »oiuw vi aq 

boom can be stabilized by torsional damping. 


interval defined by the two positive 
roots. In other words, if the torsional 
damping is gradually increased, the boom 
system will become stable, but a point will be reached when further increase 
of torsional damping will result in instability. In the right hand region 
there is always a real negative root and a pair of complex conjugate roots. 

In this region, therefore, torsional damping will never stabilize the boom 
system. 
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It follows from figure 8 that boom systems characterized by values of X Q 
close to unity can only be stabilized by torsional damping for very low values 
of A 0 * (long thermal time constants). It is clear that, in general, torsional 
damping is less effective than flexural damping. In fact, figure 7 shows that 
there is no region of the X, A* region which cannot be stabilized by flexural 
damping. 


ANALYSIS OF THE EFFECTIVENESS OF A 
PASSIVE TIP DAMPER 


One way to damp the boom system and, therefore, possibly to stabilize it 
is to incorporate a damping mechanism in the end mass. A relatively simple 
scheme is to make the end mass in the form of a thin closed shell containing a 
viscous fluid and a ball. A disturbance of the boom system then causes rela- 
tive motion between the spherical shell, attached rigidly to the boom, and the 
ball residing inside. The motion of the ball through the fluid then dissipates 
the disturbance energy. Since the ball is free to both rotate and translate, 
the damping scheme provides both torsional and flexural damping. However, it 
has been shown that torsional damping is not particularly effective in stabi- 
lizing the boom system. Therefore in this analysis it is assumed that the 
ball motion is one of translation only and the damping effectively confined to 
the flexural mode. The analysis is further simplified by assuming that 
1^, = = Dz = 0 (where, it will be recalled, and D z refer to structural 

damping only) and the mass of the damper fluid is negligible. 


The equations of motion of the damper ball along the y and z axis are 

dy d (t) 


m. 


fd 2 y(t) d 2 y d (t)“| 

l - dt* - + dt2 J 


- K 


dt 


m d 


d 2 z (t) + d2z d( t ) ~| = k dz d (t) 
dt 2 dt 2 dt 


(82) 


where 

m^ mass of the damper ball 

y^(t),zj(t) coordinates of the center of mass of the damper ball measured 
from center of mass of hollow container 


k damping constant (a force of kv is required to move the 

damper ball through the fluid with velocity v) 

The boom alone (less damper ball) force equations as given by equations (18) 
and (19) must be changed to include the forces exerted by the damper ball. 
Thus they become 


Lj 
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(83) 


Fy(t) = 


-m c 


d 2 y(t) dy d (t ) 

dt 2 dt 


F z (t) = -m c 


d 2 z(t) dZ(j(t) 

dt 2 dt 


(84) 


where m s is the mass of the damper shell fixed rigidly to the end of the 
boom. 

Combining equations (82) to (84) with equations (10), (11), (12), and 
(17) gives the characteristic polynomial which, like equation (21), has two 


factors. 

As before one of these 

is 

damped 

and need not 

be considered 

further. 

The other 

factor 

is 








m s s '{lET* 

lA ff 3 (kL) 1 -ll 

SR 2 [ 4C 3ElJJ 
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/l 3 I s 
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riiliil. 1 1) 
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m s s 2 J 

'fl(kL) 1 ] L 3 

2C EIJ 3a 

i 

0 

-k r f ^ kL) 

_ _L1 hi 

EIJ 3R 

= 0 

(85) 
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f 3 (kL) 1 I L b > 
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If P(s, m s ) is the polynomial in s formed from the first three rows 
and three columns of the above determinant, then it follows that equation (85) 
can be expressed in the form 

smurfs, m s ) + KP(s,m s + m^) =0 (86) 

It should be noted that P(s,m s ) is identical to the left hand side of equa- 
tion (42) when mj is replaced by m s and D z set equal to zero. Equa- 
tion (86) can therefore be written in the form 

smjfmsAs 3 + m s X(A - B)s 2 + s + X] 

+ K[(m s + m^)As 3 + (m s + m £ j)X(A - B)s 2 + s + X] = 0 (87) 

It is convenient at this stage to write this equation in terms of the 
dimensionless parameters given by definitions (47), (48), and (50) (where 
m T = m d + m s ) pl us those defined as follows: 


nr 



( 88 ) 


* A /m T (A - B ) 
m d 


(89) 
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Equation (87) then becomes 




m d*) „*3 


s* + (1 - mj*)A*s* + s* + X 1 


+ s* X* + s* + -A* 


'1 " m d* \ j.4 


1 - X 


) s* 4 + (1 - m d *)X* + 1 K _ x 


+ (1 + K*X*)s* + (X* + K*)s* + K*X* = 0 


Since by definition 1 - m^* 0, 1 - X >_ 0, and k* >_ 0, the coefficients of 

equation (91) are always positive. Therefore the boom system is always stable 
provided inequality (80) is satisfied. Substituting the appropriate value of 
the coefficients from equation (91) into inequality (80) yields 

3 

K * X _*Y - K* 2 {m d * + X* 2 [m d * - (1 - m d *)X]} 

- K*A*[m d *A* 2 (l - m d *) (1 - X) + m d * - (1 - m d *)X] + (1 - m d *)A* 2 X < 0 (92) 

provided X f 1 

Equation (92) may be used to answer the following question. Given a 
fixed value of X and X*, what is minimum value of m j* (subsequently denoted 
by m d m *) that will permit the boom to be stabilized? This is equivalent to 

finding the minimum value of m^* for which a real positive value of k* 
exists which satisfies equation (92) . Corresponding values of m^ * and k* 
for various values of X* and X have been calculated on a digital computer. 
The results are shown in figures 9 and 10. At md* equal to m^ * 


T T r ~n 

x X Eq(E26) 


\Boorrk 1 Boom\ 
\Boon\'c' \ D' 
\'B' \\ \ 


V/ / Values of 
\ AcandX’for 

^ ( md m) max 

I 111! I — t-+-U4 


Figure 9.- Variation of minimum m^* for stability, 
with X*, for various values of X. 


Figure 10.- Variation of k for minimum m^*, 
with X*, for various values of X. 
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equation (92) has one real negative root and a pair of equal real positive 
roots. For values of mj* below the m d m * shown in figure 9, equation (92) 

has one real negative root and a pair of conjugate complex roots. In this 
case there is no real positive value of k* that satisfies equation (92) and 
therefore the boom cannot be stabilized. It is interesting that for a fixed 
value of X, there is a value of A* for which ™d m * is a maximum. With 

this value of m^*, therefore, a boom with a given value of X can be stabi- 
lized irrespective of the material from which it is made. An analytical 
determination of the maximum value of m^*, f ° r a given X, along with the 
value of A* at which it occurs, is given in appendix E. The results are 
summarized below 


mn *) = X 

< m ' max 

(93) 

A* - 1 

(94) 

A - x 

K* = /I - X 

(95) 


It is clear that, provided the value of X associated with a given boom 
is not too large, the damping scheme provides a practical solution. However, 
booms with open cross sections can have values of X close to unity and unless 
the value of A* is either very large or very small the technique may not 
offer a practical solution. 


STABILITY OF TYPICAL OPEN SECTION BOOMS 


The analytical results given in the previous sections show that there is 
a class of booms which can exhibit severe instabilities. In order to present 
the results in a compact form they are expressed in terms of dimensionless 
parameters. This method of presentation makes it difficult to visualize the 
physical dimensions and elastic properties of booms likely to experience 
instabilities. To provide data points relative to which other booms may be 
judged and to provide convincing evidence that existing boom systems can be 
unstable, the analytical results will be applied to specific examples. The 
pertinent parameters for these examples are given in the table. These param- 
eters are typically those for silver plated copper-beryllium booms with a 
structurally open cross section. The thermal radiation constant is appropri- 
ate to a situation in which the sunline is normal to the boom axis. Both the 
radius of curvature and thermal time constant have been calculated assuming 
the boom cross section is thermally seamless (eqs. (4) and (5)). Boom A 
represents the 0G0 IV boom that is thought to have exhibited instability. 

Boom B is typical of one of the booms of a Naval Research Laboratory (NRL) 
gravity stabilized satellite which has experienced anomalous behavior. Boom C 
represents one of the booms of the Department of Defense Gravity Experiment 
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(DODGE) satellite which was observed by television camera but did not 
oscillate. Boom D represents one of the main booms of the ATS-D satellite. 

The theory given in this report is based on the assumption that the mass 
of the boom is negligible compared with the mass of the tip. This is certainly 
not true for boom A. However, if it is assumed that the mode shape has only 
a secondary influence on stability then the boom mass can be represented as an 
equivalent tip mass. The magnitude of this tip mass is such that it gives the 
correct natural frequency of the first bending mode. Thus if the boom has no 
true tip mass, as is the case for boom A, then the natural fundamental fre- 
quency is 3.52 /eTT^T where a m is the mass per unit length of the boom. 

If Amp is the equivalent tip mass then 


/ 3EI _ 3 52 /El" 

/ Am T L 3 ' /^L J 

or 

Amp - 0.242 a m L (96) 

Similar considerations in torsion lead to an equivalent tip rotational inertia 
given by 

nr = 'it rrr 

/ A^ ' 2 / OjL 2 
or 

AI^ = 0.406 ajL (97) 

Equation (97) is only an approximation, since the effect of the warping stiff- 
ness has been ignored. However, the calculation is used only to demonstrate 
that inequality (38) is satisfied (see table) and, therefore, that rotational 
inertia effects on stability can be ignored. Equation (97) is sufficiently 
accurate for this purpose. 

The values of A* and X for each of the example booms are plotted in 
figure 6 . It is clear that boom A is the most unstable, having a value of 
p* = 0.14. Therefore, in the absence of structural damping, boom A has an 
unstable oscillation which diverges as if it had a damping factor of -0.14. 1 
The value of p* for boom D is about half that for boom A. However, this 
fact does not give a balanced picture of the relative severity of the instabil- 
ities of these two booms. Thus, for example, the time to double the amplitude 
for boom A is only 0.46 minute while for boom D it is 17.20 minutes. 

This difference, of course, merely reflects the widely different frequencies 
of the two booms. However, it is important to note that in space applications 
the orientation of the run relative to the boom axis will generally vary with 
time as the satellite moves around its orbit. Therefore the conditions favor- 
able to promoting instability in a given boom may occur for only part of the 

■^When C/EI is very small /m^(A - B) - /mpL 3 /3EI therefore 
5 z f s * -P*. 
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orbit. It follows that the greater the time to double the amplitude, the less 
will be the probability that its amplitude will be severe enough to influence 
the satellite or even be detectable. This may be the reason instabilities 
were never observed in boom C. 

It follows from figure 7 that all the booms can be stabilized by flexural 
damping. In each case the required flexural damping factor is slightly higher 
than the corresponding value of |p*|. On the other hand the values of A* 
for all the booms is considerably higher than 1.0, and it follows from 
figure 8 that none of the booms can be stabilized by torsional damping alone. 

Although the booms can be stabilized by flexural damping the amount 
required for all the example booms is considerably higher than can be provided 
by the structure. The ball damper, however, seems to offer a practical solu- 
tion to this problem for booms B, C, and D. In particular, boom D could be 
stabilized with a ball damper in which the ball is about 30 percent of the 
entire tip mass (see fig. 9). In fact a reasonable design approach would seem 
to be to assume that the boom has no structural damping and to size a ball 
damper to make the boom neutrally stable. The difficult question of deciding 
the exact value of the structural damping is thereby avoided. 


CONCLUSIONS 


A linearized dynamic stability analysis of slender booms, illuminated by 
thermal radiation and subject solely to self-induced thermal bending moments, 
leads to the following conclusions. 

1. In the absence of structural damping any initially straight boom with 
either seamless or open cross section is unstable. 

2. For most practical boom systems, including all that have been flown, 
the rotational inertia about the boom axis has a negligible effect on boom 
stability. Subsequent conclusions are all based on this simplification. 

3. Initial curvature may increase or decrease the instability depending 
on its magnitude and direction. If the equilibrium shape of the boom, repre- 
senting the combined effects of initial and thermally induced curvatures, is 
such that the boom and radiation source are on the same side of a plane pass- 
ing through the boom root and perpendicular to the radiation rays, then the 
boom will be stable. In the absence of structural damping the converse is 
also true. 

4. In the absence of structural damping, booms with a seamless cross 
section are only mildly unstable. They would probably be stabilized by their 
natural structural damping. 

5. In the absence of structural damping, booms with open cross sections, 
characterized by a low ratio of torsional to bending stiffness, can be very 
unstable. Booms with an initial curvature may exhibit either oscillatory or 
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nonoscillatory instabilities depending on the magnitude and direction of the 
initial curvature. Initially straight booms, on the other hand, can only 
exhibit oscillatory instabilities. The worst case oscillatory instability has 
a frequency 66 percent of the natural flexural frequency and diverges as if it 
had a damping factor of -0.354. The magnitude of the instability is only 
slightly dependent on whether or not the cross section at the free end of the 
boom is restrained from warping. These booms can be stabilized by flexural 
damping but, in general, the amount required is greater than can be provided 
by the structure of existing booms. 

6. Torsional damping is far less effective than flexural damping in 
suppressing the instability. There exists an important class of booms of open 
cross section which cannot be stabilized by torsional damping alone. 

7. For a large class of booms, a passive ball type damper, attached to 
the tip of the booms seems to offer a practical solution to the problem of 
suppressing the instability. 


Ames Research Center 

National Aeronautics and Space Administration 

Moffett Field, Calif., 94035, November 28, 1969 



APPENDIX A 


THERMAL BENDING MOMENTS IN A THERMALLY SEAMLESS THIN- WALL BOOM 
OF CIRCULAR SECTION HEATED BY THERMAL RADIATION 


The primary aim of this appendix is to derive the linearized dynamic 
relationship between the thermal bending moment and the longitudinal twist of 
the boom. 

First the normal modes and time constants of the temperature distribution 
are determined for a boom heated by parallel rays of thermal radiation. These 
modes are then used to obtain the steady-state temperature distribution around 
the boom cross section. Knowledge of the temperature distribution then per- 
mits the steady-state bending moments acting along the length of the boom to 
be calculated. It is then assumed that the boom is suddenly rotated about its 
axis through a small angle. The accompanying variation of temperature dis- 
tribution, as expressed by the normal modes and time constants, is then used to 
determine the differential equation describing the variation of thermal bend- 
ing moment with time. To simplify the analysis, assumptions similar to those 
used in reference 1 are made. These are given below for the sake of 
completeness . 

1. The boom is sufficiently slender that heat conduction along its 
length may be neglected. 

2. The entire cross section loses heat by radiation. 

3. The radiant heat absorbed by an element on the sunlit side of the 
boom is proportional to the cosine of the angle between the surface normal at 
the element and the sun f s rays . 

4. Internal radiation is neglected. 

The first of the above assumptions permits a reduction in the dimensionality 
of the heating problem since it implies that the thermal behavior of any sec- 
tion of the boom is independent of the thermal behavior of adjacent sections. 
Thus, the entire thermal behavior of the boom can be deduced when the 
behavior of a typical longitudinal section of small length is known. 


TEMPERATURE DISTRIBUTION 


) 


» 


s 


The unsteady heat flow equations defining the temperature distribution 
around the circumference at a typical section along the boom and consistent 
with the above assumptions is 


Kh 


3 2 T (s , t ' ) 
3s 2 


pch 


3T(s ,t ' ) 
3t ' 


oeT 4 (s ,t ' ) - y (s) 


(Al) 
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along with boundary conditions representing the fact that for a seamless 
circular cross section the temperature and its first derivative are 
continuous 


T(s,t') = T(s + 2irr,t') 


(A2) 


it 


where 


3 T(s,t* ) = 3T (s + 2Trr > t’) 
3s 3s 


* 




f 


K thermal conductivity 

h cylinder thickness 

p mass density 

c specific heat 

a Stephan-Boltzmann constant 

e emissivity 

s distance measured along the circumference 

t T time 

T(s,t ! ) absolute temperature at point s at time t f 
r cross-sectional radius 

y(s) input heat flux rate due to thermal radiation 

The analysis may be simplified further by making use of the observation that 
booms are usually made from materials whose thermal properties cause the tem- 
perature difference between any two points on the boom to be small compared 
with the absolute temperature of any point on the boom, Thus, if T Q is some 
average temperature and T(s,t f ) is the difference between the temperature at 
the point s at time t r and the average temperature, then the ratio 
T(s,t f )/T 0 is usually small. This permits use of the so-called linearizing 
assumption in which squares and higher orders of T(s,t ! )/T 0 are neglected. 

By definition 


T(s,t’) = T q + T (s , t 1 ) (A3) 

and substituting for T(s,t ! ) from equation (A3), into equations (Al) and (A2) 
and using the linearizing assumption gives the following equations for 
T(s,f) 
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a 2 T(s,t*) pc 3T(s,t’) 4aeT 0 3 ^ ^ _ y(s) + 

^2 K 3t' ' Kh 1 ’ j Kh 


T(s ,t ' ) = T (s + 2irr,t') \ 


3T(s,t») _ 3T(s + 2irr,t * ) 
3s 3s 


(A4) 


(A5) 


The homogeneous part of equation (A4) is 
a 2 T(s,t') pc 3T(s ,t ' ) 


4aeT 0 3 


as" 


at' 


Kh 


T(s,t ' ) = 0 


(A6) 


A solution to equation (A6) has the form 


At, 


T(s,t' ) = e A(s) 


(A7) 


where A is a constant. Substituting for T(s,t') from equation (A7) into 
equation (A6) gives the following ordinary differential equation for the 
function A(s) 

d 2 A(s) 


ds" 


f2 2 A(s) = 0 


(AS) 


where 


4oeT 3 , 

2 _ o pcA 


= 


Kh 


(A9) 


The solution of equation (A8) is 


/.f fts fts 

A(s) = Cie + C 2 e 


(A10) 


The value of ft may now be found from the boundary conditions given by equa- 
tions (A5). Substituting for A(s) from equation (A10) into equation (A7) and 
for T(s,t ! ) from equation (A7) into equations (A5) gives 


n fts , n ^ 
Ci e + Coe 


-fts 


= Ci e 


ft (s + 27rr) 


+ C 2 e 


-ft (s+2ur) 


„ fts n — fts n ft(s + 27rr) n -ft(s + 2Trr) 
C^e - C 2 e = C^e v - C 2 e 


(All) 

(A12) 


Equations (All) and (A12) are satisfied for nontrivial values of Ci and C 2 
if and only if 


= 0 


/_ ft27ir\ 

/_ -ft27rr\ 

(l - e ) 

(i - e ) 

/, fi2irr\ 

/, -fi 2 irr\ 

(l - e ) 

-(! - e ) 


5 


1 




* 
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or 


i 


5 


s 


I 


/-. ^2ttt 

V 1 - e 

Equation (A13) has a double root e 
by 


)(i - 

ttZ-UT 


n 


e“ Q2Trr ) = 0 
1 which, in turn, 

= 0,1,2,. . . 


(A13) 

has roots given 
(A14) 


where 

i = /T 

The time constants A n can now be determined from equation (A9), 

* n = ifilai , _L (zf 

n hpc pc \r/ 

and the functions A n (s) from equation (A10), 

Ms) - C ln e in ( s/r) * C 2 n e- in t s ^ 


(A15) 


which can be rewritten in the form 

A n (s) = A ln (s) + A 2n (s) (A16) 

where 

A-m( s ) = (c ln + C2 n ^cos 

A 2 n(s) = i(c ln - C 2n )sin 


The functions Am(s) and A 2 n (s), n = 0,1,2,. . . are normal modes of the 
temperature distribution. 

It now follows that the general solution of the homogeneous equation (A6) 
is 

CO 

T(s ,t ' ) = XX nt ' [Ain (s) + A 2n (s)] (A17) 

n=0 


Equation (A17) suggests a solution to the inhomogeneous equation (A4) of the 
form 

oo 

T(s,t’) = ^ [q ln (t')cos + q 2n (t')sin (A18) 

n=0 
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where the constants Cm + ^2n an< i i(Cm - C 2 n) are included in the 
functions q ln (t') and q 2n (t'), respectively. Substituting equation (A18) 
into equation (A4) and using equations (A8) and (A9) yields the equation 


E lpqin(t') 




ns 

cos — + 
r 


d v (t,) 

tit ' + j 


YCsj-aeTp 4 

hpc 

(A19) 


which must be satisfied by q ln (t') and q 2n (t') for equation (A18) to be a 
solution. 


i 


f 


STEADY- STATE BENDING MOMENTS 


Consider now the situation as shown in figure 
axes are assumed to be fixed in the cross section. 


11(a), where the y M and z" 
By assumption 3 the input 


Thermal radiation 


s 



Thermal radiation 



Figure 11.- Coordinate systems for thermal analysis. 


heat flux rate due to solar radiation is given by the expression 

Y (s) = J s ot s 0(s)cos ^ 


(A20) 


where 

J s = thermal radiation constant 
a s = absorptivity 

and 

0 (s) = 1 for cos — > 0 

r “ 

= 0 for cos — < 0 
r 
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The variation of y(s) with s/r is 
shown in figure 12. It follows from 
equation (A20) that y(s) is an even 
function of s/r. Therefore, since 
aeT Q 4 is a constant, the left hand side 
of equation (A19) is also an even func- 
^ tion of s/r. Since the left hand side 
of equation (A19) is a Fourier series 
with argument s/r, it follows that equa- 
Figure 12.- Circumferential distribution of input tion (A19) can be satisfied if and only 
heat flux rate. if q 2n (t f ) = 0 for all n. 

It is shown in reference 1 (after sign adjustments to account for the 
adopted axis system) that the thermally induced bending moments about the 
y M and z M axes are given by the expressions 


Myii ft ' ) = e c Eh J 

J 27rr T(s,t')z"(s)ds 

(A21) 

= -e c Eh 

/ 2lTr T(s s t ' )y" (s)ds 

(A22) 


where 



e c coefficient of thermal expansion 

E Young f s modulus 


y"(s)l 

z"(s)J 


y M and z" coordinates of a point s on the perimeter of the boom 


It follows from figure 11 that 


y"(s) = 

s 

-r cos — 
r 

(A23) 

z"(s) = 

s 

-r sm — 
r 

(A24) 


Substituting equation (A18) into equations (A21) and (A22) (noting that 
q 2n (‘ t? ) E 0 for all n) and performing the integrations indicated gives the 

following expressions 


Ct') = 0 

(A25) 

(f) = e c EhTrr 2 qil (t') 

(A26) 
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However, it follows from equation (A19) that [dq^ (t ' ) /dt * ] + X^q^Ct') is 

equal to the coefficient of cos (s/r) in the Fourier expansion of 
y(s) - oeT 0 4 /hpc. Thus, 


dq n (t') 

dt' 


+ 


* iq il 


J s a s 

2hpc 


(A27) 


Equations (A26) and (A27) may now be combined to give the following expression 
for M „(t') in terms of the material properties. Thus, 


dM », (t ' ) 
dt 1 + 


MM z „(t') 


J s a s e c E^r 

2pc 


2 


(A28) 


Equations (A25) and (A28) show that in the steady state, 


My.,, (“0 

M Z „C») 


= 0 

J s a s e c E^r 2 

2pcXj 


(A29) 




( 


CHANGE OF BENDING MOMENTS DUE TO SMALL ROTATION 


Suppose, now, that the boom, whose cross section is shown in figure 11, has 
reached its steady-state temperature distribution. The bending moments acting 
on the boom are given by equation (A29) . Suppose, further, that the cross 
section is instantaneously rotated through a small angle ip . It is required 
to know the equations defining the subsequent change in the thermal bending 
moments . 


The equation governing the perturbation in temperature caused by the 
small rotation differs from equation (A4) only in that the right hand side is 
replaced by -Ay(s)/Kh where Ay(s) is the change of heat flux rate. It should 
be noted here that t is the time measured from the instant of rotation 
through the angle ip. The equation corresponding to equation (A19) is, 
therefore 


oo 

E 

n=0 


_ dSq m (t) 

at * VWjnCU 


ns 

cos — + 
r 


dAq 2n 


dt +A n A q 2n (:t:) 


. ns I _ Ay(s) 


sm 


hpc 


(A30) 


The instantaneous change of heat flux rate is given by (see eq. (A20)) 


where 


Ay(s) = J s a s [e 1 (s)cos (7 + - © 2 (s)cos 

Qj (s) = 1 C0S (7 + ^ ^ 0 


(A31) 


= 0 cos (f' + ^ < 0 




t 
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0 2 (s) = 1 cos — >_ 0 

= 0 cos — < 0 

r 

If <p is assumed to be sufficiently small that only terms linear in ij > need 
be retained, equation (A31) becomes 


Ay (s) = -J s a s 0(s)ijj sin — 


(A32) 


The variation of Ay(s) with s/r is 
shown in figure 13. Clearly, Ay(s) is 
an odd function of s/r; therefore, equa- 
tion (A30) shows that Aq ln (t) = 0 for 
all n. 


\. Proceeding as for equations (A25) 

N and (A26) gives corresponding expres- 

Figure 13.- Change of circumferential distribution sionS for the changes in the thermal 
of input heat flux rate due to twist. bending moments; thus. 


AM „(t) = -e c Ehnr 2 Aq 21 (t) 


(A33) 


AM z „(t) = 0 


Equation (A30) shows that the expression [dAq 2 1 (t) /dt] + A 1 Aq 2 (t) is equal 
to the coefficient of sin(s/r) in the Fourier expansion of Ay(s)/hpc. Thus, 


dAq (t) 
dt + 


J s a s^ 

2hpc 


(A34) 


Equations (A33) and (A34) may now be combined to give the following expressions 
for AM (t) 


dAM ,, (t) 


dt + A l AM y-'( t ) = 


J s a s e c Eirr 2 i|) 

2pc 


(A35) 


Equation (A35) may now be combined with equation (A29) to yield 


dAM (t) 

Jt + X l AM y u Ct) = X l M z" ( °°^ 


(A36) 
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AM ,,(t) may now be resolved about a set of axes y* and z ' 
11(b)) which are fixed relative to the direction of the incident 


The bending moment 

(see fig 
radiation. The component 


y' axis is given by 


AMy, (t) about the 

AM ,,(t)cos ip ~ AMy, (t) while the component AM z ,(t) about the z' axis is 
given by AM^,,(t)sin ip ~ 0 (since by eq. (A36) AM^, (t) is the same order as 


\p) . The overall conclusion, therefore, is that following an instantaneous 
rotation ip, the thermal bending moments acting on the cross section consist 
of 


(a) A moment M^,, (°°) whose direction rotates with the cross section (body 
fixed moment) and whose magnitude is constant and given by equation (A29) 

(b) A moment AMy, (t) whose direction is along the rays of the incident 
radiation and whose magnitude is given by 

dAM , (t) 

^ + A!AM yI (t) = A 1 M z , I (oo)^(t) (A3 7) 

where iKt) is regarded as a function of time. 
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APPENDIX B 


DEFLECTION EQUATIONS FOR A CIRCULARLY BENT BOOM 
OF CIRCULAR CROSS SECTION 



b(L) 


Figure 14.- Coordinate system for deflection 
analysis . 


Consider a boom whose locus of flexural centers, in the unloaded condi- 
tion, forms the arc of a circle of radius R. Furthermore, let the boom be 
attached rigidly at one end to form a cantilever. It is assumed from the out- 
set that the length of the boom L is small compared with the radius of 
curvature R. The situation is shown in figure_14 which also shows the 
coordinate systems adopted. The unit vectors t(x), iT(x) , and b(x) at a 

typical point A are the tangent, 
normal, and binormal, respectively, to 
the locus of flexural centers. The 
coordinates of point A relative to the 
t(0), n(0), b(0) axes are x Q , y Q , z Q 
while the deflection_ of point A rela- 
tive to the n(x) , b (x) axes are denoted 
by y(x) and z(x). The angle of twist 
about t"(x) is denoted by iKx) . 

The problem is to determine the tip 
deflections y(L) and z(L) and the twist 
4»(L) due to the simultaneous action of 
forces Fy and F z acting along n(L) 
and b(L), a couple T acting along 
t(L) and an arbitrary distribution of 
bending moments AM n (x) and AM^Cx) 
acting along n(x) and b(x). The loads 
on the boom are shown in figure 15. The 
assumption is now made that the boom 
deflections are sufficiently small that 
the principal of superposition holds and 
that, therefore, the bending moments and 
torques remain constant and equal to 
their values calculated as if the boom 
were undeflected. It follows from fig- 
ure 15 that the total bending moments 
M n (x) and (x] along n(x) and b (x) , 
respectively, at point A, are 

M n (x) = -F Z BC + T sin (^X) + AM n (x) 

M b( x ) = FyAD + AMb(x) 

while the torque M t (x) along t (x) is 
given by 



Figure 15.- Loads acting on the boom. 
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The distances 


M t (x) = F Z AC + T cos — ^ — 

BC, AD, and AC are given by the following expressions 
BCf = AD = R sin L n X 


AC = R 




COS 



X 


) 


It is assumed at this point that the order of magnitude of the ratio of the 
couple T to either force F z or Fy never exceeds L 2 /R. Expanding the 
above expressions into powers of (L - x)/R and neglecting all above the third 
power in the moments due to Fy and F z and the first power in moments due to 
T (in agreement with the above assumption) yields the following expressions 
for M n (x) , M b (x), and M t (x) 


M n 00 

M b (x) 

M t (x) 


-F 


?x 


z 


|^(L - x) 

(L - x) - 
(L - x) 2 


2R 


- - L 6 ^ 2 X ' ) j + ^ (L - x) + AM n (x) 

( - L 6 ^ 2 X - ) j + AMb(x) 

+ T 


(Bl) 

(B2) 

(B3) 


It is convenient to start with the calculation of tJj(x) . The twist angle 
of a bent boom can be regarded as the sum of two components, one representing 
the contribution due to twisting torque, the other the contribution due to 
bending. Thus, 


iK*) = <l> T (x) + ^ B 00 (B4) 

The relationship between angle of twist and applied torque is given by the 
differential equation (ref. 6) 


“Ci 


d 3 t|> T (x) 

dx 3 


# T (x) 


c T = M tW 


(B5) 


The appropriate solution of equation (B5) must satisfy the following end 
conditions 


<M X ) 


x=o 


= 0 


(B6) 
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x =0 


0 


(B7) 


d^ T (x) 


dx 


along with 


d^OO 


or 


dx 


d 2 ^ T (x) 


= 0 i£ Tip Warping is Restrained (TWR) 


x=L 


dx 2 


= 0 if Tip Warping is Unrestrained (TWU) 


x=L 


The appropriate solution of equation (B5) is given by 


where 


/:(• 


<Kp(x) = I jB sinh(kx) - f M t (x)sinh[k(x - x)]dxldx 

1 •'o I 




g = — — — f M t (x)sinh[k(L - x)]dx (TWR) 

C 2 k sinh(kL) Jq 


or 


g = 


Cjk cosh(kL) J q 




(x)cosh[k(L - x)]dx (TWU) 


Substituting equation (B3) into equations (BIO), (Bll), and (B12) and 
performing the integrations gives the following results 

t T W = ^ (* 3 + 3Lx(L - x) + | | C s^h( + kL ~ j ' {cosh f k (L- - cosh(kL)} 


A 6[1 - cosh(kx)] \ T 
k 3 sinh(kL) ) + C 


cosh(kx) - cosh[k(L- x)] 
k sinh (kL) 


tanh 


CBS) 

(B9) 

(BIO) 


(Bll) 

(B12) 


(TWR) 

(B13) 
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or 


♦ T oo 


- — (■ 

6 RC \ 


x d + 3Lx(L 


x)+ ^l + iimi±2i {sinh[k(L . x)] 

k 2 k 3 cosh (kL) 


sinh (kL) } 


) 


T ( sinh[k(L-x)] tanh(kL) 
+ C j x + "k 'cosh(kL) ‘ k 


(TWU) 


When x = L 


where 


or 


. , n _ F z L 3 fi(kL) TLf 2 (kL) 

n*- LJ “ 6 RC C 


f 1 (kL) = 1 + 


£ 2 (kL) = 1 


6 _ 3 [ (kL ) 2 + 4] 

(kL ) 2 (kL ) 3 


tanh 


(t)1 


TWR 


kL 


tanh 




£l (kL) = 


f 2 (kL) = 


. 6 3[(kL ) 2 + 21 , ,, . , 

1 + = — - — 5 — - tanh(kL) 

(kL ) 2 (kL ) 3 


1 - tanh (kL) 


TWU 


(B14) 


(B15) 


(B16) 


Graphs of fi(kL) and f 2 (kL) are given in figure 3. 

The twist due to bending ^g(x) may be obtained from the strain energy 
due to bending by application of Castigliano ' s theorem (ref. 7). It is first 
necessary to introduce, in addition to the loads indicated in figure 15, a 
dummy torque T t acting at the point x and along the tangent to the boom 
axis. The bending moment at the point x about the n(x) axis then becomes 


M n '(x) = M n (x) + T t sin — ) 

= Mn(x) 


X >_ X 


X < x j 


(B17) 


while the bending moment about the b (x) axis remains unchanged. The strain 
energy in bending is given by the expression 


U B 



[Mn ' (x) ] 2 [Mb (x) ] : 


2EI 


2EI 


dx 


(bis) 
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and by the Castigliano theorem 




M n(x) 3 [M n * (x)] 


El 


3T+ 


dx 


(B19) 


Differentiating equation (B17) with respect to T t and substituting for 
3 [M n ' (x) ] /3Tt in equation (B19) yields 


>I> B 00 = ¥T / M n (x)sin ( X ~ X ) dx 


(B20) 


Substituting for Mn(x) from equation (Bl) into equation (B20) and performing 
the integration indicated yields the following expression for t(< B (x) 


^r(x) 


EIR 


(^T _ ?pj + J AM n (x)sin(^-^) dx (B21) 


where the simplifications permitted by the assumption made about the relative 
magnitudes of the forces Fy and F z to the torque T have been made and 
terms of order F Z L 5 /R 3 and higher have been neglected. 

Combining equations (B15) and (B21) yields 

- A] * * A J o L (HA) * 

(B22) 

The linear deflections y(L) and z(L) may also be regarded as the sum of 
two components. In this case they are the contributions due to the bending 
moment and to the angle of twist ^p(x) 


y(L) = y B (L) + y T (L) 
z (L) = Zg(L) + z t (L) 


(B23) 


The contributions y B (L) and z B (L) may be calculated from the strain energy in 
bending by an application of Castigliano f s theorem. Thus the stain energy in 
bending is given by equation (B18) where in this case the dummy torque T t 
need not be included. The expressions for y g (L) and z B (L), by Castigliano, 
are 




CL) 


3U B 


3F, 


1 

El 



3Mb (x) 
3F„ 


dx 


(B24) 
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dx 


(B25) 


z b (L) 


3U b 

3F z 


1 

El 



M n (x) 


9M n (x) 

3Pz 


where use has been made of the facts that, from equations (Bl) and (B2) , 
9M n (x)/3Fy = 3M^(x)/9F z = 0. Substituting from equations (Bl) and (B2) into 
equations (B24) and (B25) yields 




[ 


1 ^ 


] dx + w/ AM b (x)sin (^) 


dx 


(B26) 


z b C l ) 


-£J l [' 


(L- x) 


(L - x) 3 
6R 2 


■ dx - J (L - x) £(L - x) - 


(L - x) 3 
6R 2 


dx 


R 

El 


J AM n (x)sin ^- L - dx 


(B27) 


where the quantity (L - x) - [ (L - x) 3 /6R 2 ], which is an approximation to 
R sin[(L - x)/R], has been replaced by its true value in the two integral 
expressions. Equations (B26) and (B27) simplify to 


y B ( L ) 

z B (L) 


r y 

3EI 

F z 

3EI 


( l3 ‘ IP - ) * ET { “b Cx)sin (^-p) dx 

( l3 - 1^2-) ■ Sir-CT J f 4 M nW sin ( i i ii ) dx 


(B28) 

(B29) 


It follows from the geometry of the boom in the deflected state that, to 
the required accuracy, the deflection y T (L) , due to twist, is zero, while the 
deflection z T (L) is given by the expression 



Figure 16.- Geometry of the linear deflection due 
to twist. 


(see fig. 16) (B30) 

Integrating equation (B30) by parts 
gives 

z T (L)=jf ^ T (x)sin dx (B31) 

Substituting for ipj(x) from equa- 
tions (B15) and (B16) into equation (B31) 
and approximating sin[(L-x)/R] by 
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(L - x)/R yields 


F z L 5 f 3 CkL) TL 3 f! (kL) 
20R 2 C 6RC 


where 


,, T . , 40 20 5 

3( ) + 3(kL) 2 + (kL) 4 ~ (kL) 5 


J 4 [ (kL) 2 + 2]tanh^y) + (kL) 4 coth(kL)J TWR 


40 20 5 [ (kL) 2 + 2] 2 . . 

f3(kL) = 1 + i(n^ + ' (kT ) 3 ~ tanh(kL) 


Graphs of f 3 (kL) are given in figure 3. 

It follows from equations (B28) , (B29) , and (B32) that 


y /r 3 _ l: 


z f ™ b ( x > sin (t) 


, f Li - P riiliw 

ztL - ) F z 1 3EI 5R 2 [ 4C 


~ f 3 (kL) Jj) TL^_ [ 

4C 3EI_|j 3R L" 


TL 3 f f 1 (kL) 1 "I R 


3R I 2C El El 


W I AM n (x)sin( 


APPENDIX C 


ORDER OF MAGNITUDE OF ANGLE OF TWIST DUE TO 
THERMAL BENDING MOMENTS 


At any point x along the boom, the angle of twist due to thermal 
bending moments is given by the expression [N(x,t)/RR s ]sin a, where 

N (x, t) = R f n(x,t)sin 

•'O 

These two expressions follow from equations (12) and (13) with L replaced by 
x. The quantity n(x,t) is given by equation (6), which has the general 
solution 


(H 1 ) 


dx 


(Cl) 


n(x,t) = e Xt n(x,0) + A e ^ ^ T h(x,x)dx 


(C2) 


Since, by definition, the boom is undisturbed at t = 0 it follows that 
n(x,0) = 0. Also \p(x,r) is continuous in x and is positive and 

integrable over the region 0 <. x <_ t. Therefore the first integral theorem 
of mean value may be applied to equation (C2) giving 


n(x,t) 


\ip(x,T) 



e -A(t-x) 


dx 


(C3) 


where 0 < x <_ t. Therefore 


n(x,t) = i|)(x,x)(l - e Xt ) 


(C4) 


Since n(x,t) is continuous and (x - x) is positive and integrable over 
the interval 0 <_ x <_ x, the first integral of mean value may be applied to 
equation (Cl) giving 


N(x,t) 


Rn (x,t) 


f •*» (H 1 ) 


where 0 <. x <_ x. Therefore, 


N(x,t) 


Rn(x,t) (l - cos 


(C5) 
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( 


It follows from equations (C4) and (C5) that 

N(x,t) = R<J>(x,7)(l - cos - e Xt ) 

or 

N(x,t) < R^(x,t)^ 1 - cos ^ (t > 0) (C6) 

If it is assumed that \Jj(x,7) > ^(x,7), inequality (C6) can be expressed in 
the form 


N(x,t) 

RR S 


sin a 


< 'Hx,t) 


x 2 sin a 
2RR S 


(t > 0) 


(C7) 


where the factor [1 - cos(x/R)] has been replaced by its approximate value 
x 2 /2R 2 . Inequality (C7) indicates that over any interval of time starting at 
zero the maximum angle of twist due to thermal bending moments alone is less 
than the maximum total angle of twist multiplied by the factor (x 2 /2RR s )sin a. 
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APPENDIX D 


EVALUATION OF THE REAL PART OF THE ROOT OF THE 
STABILITY POLYNOMIAL CORRESPONDING TO 
THE MOST RAPID DIVERGENCE 


The stability polynomial for the special case of 1^ = D^ = 0 is 
(eq. (51)) 



where the variables are defined by equations (47) to (50). In this appendix 
the real part of the root of equation (Dl) corresponding to the most rapid 
divergence is sought under conditions of (a) varying X* at a fixed X, and 
(b) varying X at a fixed A*. 

It follows from the stability criteria for a cubic given in reference 5 
that the unstable mode, if it exists, must be oscillatory. It is therefore 
assumed that, in the region of interest, equation (Dl) has one real root and a 
pair of conjugate imaginary roots. Thus the characteristic polynomial must 


have the following representation 

(s* + a*) (s* 2 - 2p*s* + p* 2 + q* 2 ) = 0 (D2) 

where p* is the real part of the most rapid divergence and a* and q* are 
both real quantities. Equation (D2) may be expanded to give 

s* 3 + s* 2 (-2p* + a*) + s(p* 2 + q* 2 - 2p*a*) + a*(p* 2 + q* 2 ) = 0 (D3) 

Equating coefficients of equations (Dl) and (D3) 

-2p* + a* = 2? z + A* (1 - X) (D4) 

p* 2 + q* 2 - 2p*a* = (1 - X)(l + 2A*<; Z ) (D5) 

a* (p* 2 + q* 2 ) = A* (1 - X) (D6) 


VARYING A*, FIXED X 

At a stationary value of p* with respect to X* it is necessary that 
9p*/3A* = 0. The result of differentiating equations (D4) , (D5) , and (D6) 
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with respect to A* is 

3a * - 1 - X 
3 A* 1 X 

(D7) 


i$r - 2p* ll. = 2tzCl - x) 

(D8) 

3 a* 
3A* 

(p* 2 » q* 2 ) + 2a*q* |3l = (1 - X) 

(D9) 


Equations (D7) , (D8) , and (D9) hold for nontrivial values of 3a*/3A* and 
q*(3q*/3A*) if and only if 


1 0 -Cl - X) 


2p* 

2 

-2C Z (1 - X) 

+ q* 2 

2a* 

-a - x) 


(DIO) 


If x / 1, equation (DIO) simplifies to 

p* 2 + q* 2 + 2a*p* + 2a*C z = 1 (Dll) 

It also follows from equations (D4) and (D6), after eliminating A*(l - X), 
that 

a* (p* 2 + q* 2 ) = -2p* + a* - 2z z (D12) 

p p 

and from equations (Dll) and (D12) by eliminating p* + q* that 

a* 2 = 1 (D13) 

Since p* and q* are real by definition, equation (D6) shows that a* must 
be positive for X < 1. Therefore 


a* = 1 


(D14) 


Equations (D5) , (D6), and (D14) may be combined to eliminate A*, 
yielding 


p* 2 + q* 2 - 2p* = 2? z (p* 2 + q* 2 ) + 1 - X (D15) 

2 2 

Eliminating p* + q* between equations (Dll) and (D15) gives the following 
expression for p* 


p* = 


4(1 - C Z ) 


(D16) 
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Substituting for p* from equation (D16) into equation (D4) produces the 
following equation for X* 


The imaginary part 

thus , 


1 - X [ X " 2(1 - c z ) ] 

of the complex root may be obtained from equation 
q* 2 = 1 - p* 2 - 2 (p* + ? z ) 

VARYING X, FIXED X* 


(D17) 

(Dll); 

(D18> 


At a stationary value of p* with respect to X it is necessary that 

3p*/9X = 0. The result of differentiating equations (D4) , (D5) , and (D6) 
with respect to X is 


3 a* 
3X 


= -X* 


2q*^ 


8X 


2 P* H 1 ' -n * 


9a" 

8X 


(p* 2 + q* 2 ) + 2a*q* |3l = -X* 


(D19) 

(D20) 

(D21) 


Equations (D19), (D20) , and (D21) hold for nontrivial values of 3a*/3X and 
q*(3q*/3X) if and only if 


-2p* 


p* + q* 


0 

2 

2 a* 


1 + 2X*£, 


= 0 


(D22) 


Equation (D22) simplifies to the following expression 


2p*a* = 


-a*(l + 2X*? Z ) + X’ 


(D23) 


It follows from equations (D5) and (D6) by eliminating (1 - X) that 
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(D24) 


(p* 2 + q* 2 ) [A* - a* (1 + 2A*£ Z )] - 2p*a*A* = 0 

Equations (D12) , (D23) , and (D24) are linear in (p* 2 + q* 2 ) and 2p*. The 
corresponding eliminant is 


a* 

1 


1 

a* 


-a* + 2t, z 

a* (1 + 2A*C Z ) - A* 
A* 


A* - a* (1 + 2A*£ z ) -a*A* 


0 


(D25) 


which simplifies to 

(-a* + A*) 2 = a* 2 (A* 2 - 2A*? z ) 

It follows from equation (D26) that this solution can hold only if 


A* 2 - 2A*C Z >. 0 


(D26) 


or 



The appropriate solution of equation (D26) is 


(D27) 


-a* + A* = a* A* 2 - 2A*C Z 


or 


a* = y (D28) 

1 + A* - 2A*C Z 

2 2 

Eliminating p* + q* between equations (D12) and (D24) and substituting for 
a* from equation (D28) yields the following expression for p* 


p* = 


2(1 + 


2 A~* 


- e. 


- 2A*t - 2A*£ Z ) 


(D29) 


Equation (D4) then provides the value of X corresponding to the fastest 
divergence , thus , 


X = 1 


A* 2 - 2A*£ z - 2A*£ Z 

(1 + A * 2 - 2A*? Z )(1 + 2A* 2 - 2A*^ Z - 2A*? Z ) 


(D30) 
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Since X < 1 it follows from equation (D30) that 


A * 2 - 2 A*£, - > 0 

*-* Zr 


or 


„ A + 4A* 2 - 1 
4A* 


(D31) 


Since (/l + 4A* 2 - 1)/4X* < X*/2 it follows that inequality (D31) is a more 
stringent limitation on permissible value of r >z than is inequality (D27) . 

The imaginary part of the complex root may be obtained from equation (D12) ; 

thus. 


q 


(p* 


Sz) + 1 


(D32) 
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APPENDIX E 


EVALUATION OF THE MAXIMUM VALUE OF m^* 

For any value of X and A* the value of m dm * given in figure 9 
satisfies equation (92) with m d * replaced by m d m ** Thus, 

T^r- - “* 2{,n v * X * 2 K,* - (1 - 

-K*A*[m dm *A* 2 (l-m dm *)(l-X)+m dm *- (1 - m dm *)X] + (1 - m dm *)A* 2 X = 0 (El) 

Furthermore the computer study which resulted in figure 9 shows that equa- 
tion (El) has one real negative root and a pair of equal positive roots. 


Therefore equation (El) may be represented in the form 

(k* + a) (k* - b) 2 = 0 (E2) 

or 

k * 3 + K* 2 (a - 2b) + K*(b 2 - 2ab) + ab 2 = 0 (E3) 

Equating coefficients of equations (El) and (E3) produces 

a - 2b = -{m^* + A* 2 ^^* - (1 - m dm *)X]} - 1 x - ~- (E4) 

b 2 - 2ab = -[m dm *(l-m dm *)(l- X)A* 2 + mdm *- (1 - m dm *)X) (E5) 

ab 2 = (1 - m dm *)(l - X) A* (E6) 


A necessary condition for m d m * to a maximum with respect to variations in 
A* is (9m dm */3A*) = 0. Differentiating equations (E4) , (E5) , and (E6) with 
respect to A* and using this condition gives the following equations 


(b - a) 


9a 
3 A* 

23b 

3A* 

1 - X 

2 

X* X 

{m d m * 

9b 

b 3a 

1 - 

X 

ax* 

II 

* 

1 

X 

~ m d m 


hn 


J m 


-4n 


(E7) 

(E8) 
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b 2 


9A* 


+ 2ab 


1^. - (1 - „ dm *) (1 - X) 


(E9) 


Eliminating 9a/9A* and 9b/9A* gives 


tad/ - »* 2 (»d m * - a - "U/Wb 2 


+ 2 


(1- X) 2 


m din *(l-m dm *)A*b- Cl-m dm *)(l-X) = 0 (E10) 




m 


It is convenient at this stage to introduce a variable £ defined by 

5 - ”V -x— 

Equations (E4), (E5) , (E6) , and (ElO) can then be rewritten as 

n 

1 + X* 


a - 2b = -5 


+ (1 - m dm *)(l - X) X* 


b 2 - 2ab = -?[X* 2 (1 - md m *)(l - X) + 1] + (1 - m^*) (1 - X) 


ab 2 = (1 - m dm *)(l - X)X* 

2 . 

- ~ x v* + Cl - m dm *) a - X)x*jb 2 - 2?X* 2 b(l - m dm *)(l - X) 


CE11) 
(E12) 
CE 13) 


+ (1 - m dm *) (I - X)A* = 0 (E14) 


Substituting ab 2 for (1 - mj *) (1 - X)A* (eq. (E13)) in equations (Ell), 
(E12) , and (E14) gives 111 


a - 2b = -5 - y* - — + ab 2 


b 2 - 2ab = -£(A* ab 2 + 1) + 


ab 2 


CE15) 

(E16) 


5 - '„ A * + ab 2 - 2?A*ab + a = 0 


CE17) 


Substituting E, from equation (E15) into equations (E16) and (E17) gives, 
after some simplification. 
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A* 3 (a 2 b - 2ab 2 - a 2 b 3 + 2a- b)b + A* 2 (a - 2b) + A*(2ab - b 2 ) + ab 2 = 0 (E18) 


A* 2 (a 2 b - 2ab 2 - a 2 b 3 + a - b) + (b + ab 2 ) = 0 


(E19) 


If a solution exists then equations (E18) and (E19) must have a common root. 
The dialytic eliminant of equations (E18) and (E19) is 


(a 2 b - 2ab 2 - a 2 b 3 + 2a - b)b 


a - 2b 


(a 2 b - 2ab 2 - a 2 b 3 + 2a - b)b 


a 2 b - 2ab 2 - a 2 b 3 + a - b 


2ab - b 2 


a - 2b 


a 2 b - 2ab 2 - a 2 b 3 + a - b 


a 2 b - 2ab 2 - a 2 b 3 + a - b 


b + ab 2 


ab 2 0 

2ab - b 2 ab 2 

0 b + ab 2 

b + ab 2 0 


= 0 


(E20) 

It can be shown by direct substitution that equation (E20) is satisfied when 

a = b (E21) 

Substituting a for b in equation (E18) then yields 


a 


(E22) 


From equation (E 15) substituting for X* and b from equations (E22) and (E21) , 
respectively, gives the following expression for g 


5 = m d. 


1 - X 


■m X 


Eliminating a from equations (E23) and (E13) then gives 


m . * = 

d m 


Om/) 


max 


from which it follows that 


a = A - X 

From equations (E22) , (E24) , and (E25) 


(m d *) =1 V 

v u m max \* z 


(E23) 


(E24) 


(E25) 


(E26) 


The graph of equation (E26) is shown by the dashed line in figure 9. 
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TABLE 1. EXAMPLE BOOMS 


Parameter 

Radius of cross section, 
r, m 

Wall thickness, h, m 
Thermal conductivity, K, 
J/sec m °K 

Mass density, p, kg/m 3 
Specific heat, c, J/kg °K 
Coefficient of thermal 
expansion, e c , m/m °K 
Stephan-Boltzmann constant, 
a, J/sec m 2 °K 4 
Thermal radiation constant, 
J s , J/sec m 2 
Emissivity , e 
Absorptivity, a s 
Young's modulus, E, n/m 2 
Second moment of area, I, 
m 4 

Torsional stiffness, C, nm 2 
Warping stiffness, Ci, nm 4 
Absolute ambient temperature, 
Tq > °K 

A, 1/sec (eq. (5)) 

R s > m ( e q- C 8 )) 

Mass per unit length, a m , 
kg/m 

Polar inertia per unit 
length, aj, kg m 
Boom length, L, m 
True tip mass, my, kg 
Effective tip mass, m T +Am T , 
kg 

True rotational inertia of 
tip mass, 1^, kg m 2 
Effective rotational inertia 
of tip, 1^ + AIjj,, kg m 2 

L/C/Ci = kL 


Boom A 


Boom B 


Boom C 


0.634xl0" 2 

0.5xl0“ 4 

0.126xl0 3 

0.75xl0 4 

4.18xl0 2 

0.187xl0 -4 

5.71xl0- 8 

1.385x10 s 

0.035 

0.13 

1.31xlO n 

0.41xlO" 10 

0.93xl0“ 4 

0.262xl0' 2 

300 

1.0 

600 

0.019 

0.76xl0" 6 


Boom D 


18.3 

18.3 

24.4 

37.8 

0 

1.61 

3.76 

3.63 

0.084 

1.694 

3.872 

3.804 

0 

1 . 605x10" 3 

6x10” 3 

6xl0" 3 

0.565X10" 5 

1. 611X10- 3 

6 . 008x10" 3 

6.012x10" 3 

3.44 

3.44 

4.60 

7.11 
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TABLE 1. EXAMPLE BOOMS - Concluded 


Tip condition 

Unrestrained 

Restrained 

Restrained 

Restrained 

flCkL) (fig- 3) 

0.49 

0.41 

0.52 

0.665 

£ 3 (kL) (fig. 3) 

0.29 

0.275 

0.37 

0.51 

A (eq. (27)) 

1.2691xl0 3 

1. 2232xl0 3 

5.6802x10 s 

6.2127xl0 4 

B (eq. (41)) 

8.8883xl0 2 

8.4286xl0 2 

4.7788x10 s 

5.8777xl0 4 

D (eq. (29)) 

8.9688x10 s 

7.5044x10 s 

2.2561xl0 4 

1.0725x10 s 

(1//I0)(A/D) (eq. (33)) 

0.045 

0.0515 

0.0796 

0.184 

Radius of gyration of 





equivalent tip mass. 





P ip > m 

0.00632 

0.0308 

0.0402 

0.396 

/mpfA - B) and X* 

5.655 

25.383 

58.217 

111.887 

X 

0.700 

0.689 

0.8413 

0.946 

p* (fig. 6) . 

0.142 

0.042 

0.048 

0.075 

Time to double amplitude 





(no damping) , min 

0.46 

6.98 

14.02 

17.20 
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